Densities and Distributions 



Back in section 12.1 I presented a careful and full definition of the word "function." This is 
useful even though you should already have a pretty good idea of what the word means. If you haven't 
read that section, now would be a good time. The reason to review it is that this definition doesn't 
handle all the cases that naturally occur. This will lead to the idea of a "generalized function." 

There are (at least) two approaches to this subject. One that relates it to the ideas of functionals 
as you saw them in the calculus of variations, and one that is more intuitive and is good enough for 
most purposes. The latter appears in section 17.5, and if you want to jump there first, I can't stop you. 

17.1 Density 

What is density? If the answer is "mass per unit volume" then what does that mean? It clearly doesn't 
mean what it says, because you aren't required* to use a cubic meter. 

It's a derivative. Pick a volume AV and find the mass in that volume to be Am. The average 
volume-mass-density in that volume is Am/ AV. If the volume is the room that you're sitting in, the 
mass includes you and the air and everything else in the room. Just as in defining the concept of 
velocity (instantaneous velocity), you have to take a limit. Here the limit is 

Am dm 

lim — Y7 = -TTT (I'-l) 

AV^o AV dV ^ ' 



Even this isn't quite right, because the volume could as easily shrink to zero by approaching a line, and 
that's not what you want. It has to shrink to a point, but the standard notation doesn't let me say 
that without introducing more symbols than I want. 

Of course there are other densities. If you want to talk about paper or sheet metal you may 
find area-mass-density to be more useful, replacing the volume AV by an area AA. Maybe even linear 
mass density if you are describing wires, so that the denominator is Ai. And why is the numerator a 
mass? Maybe you are describing volume-charge-density or even population density (people per area). 
This last would appear in mathematical notation as dN/dA. 

This last example manifests a subtlety in all of these definitions. In the real world, you can't take 
the limit as AA — > 0. When you count the number of people in an area you can't very well let the 
area shrink to zero. When you describe mass, remember that the world is made of atoms. If you let 
the volume shrink too much you'll either be between or inside the atoms. Maybe you will hit a nucleus; 
maybe not. This sort of problem means that you have to stop short of the mathematical limit and let 
the volume shrink to some size that still contains many atoms, but that is small enough so the quotient 
Am/ AV isn't significantly affected by further changing AV^. Fortunately, this fine point seldom gets 
in the way, and if it does, you'll know it fast. I'll ignore it. If you're bothered by it remember that 
you are accustomed to doing the same thing when you approximate a sum by an integral. The world is 
made of atoms, and any common computation about a physical system will really involve a sum over all 
the atoms in the system (e.^. find the center of mass). You never do this, preferring to do an integral 
instead even though this is an approximation to the sum over atoms. 

If you know the density — when the word is used unqualified it commonly means volume-mass- 
density — you find mass by integration over the volume you have specified. 

m= pdV (17.2) 
Jv 

* even by the panjandrums of the Systeme International d'Unites 
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You can even think of this as a new kind of function m{V): input a specification for a volume of space; 
output a mass. That's really what density provides, a prescription to go from a volume specification to 
the amount of mass within that volume. 

For the moment, I'll restrict the subject to linear mass density, and so that you simply need the 
coordinate along a straight line. 



A(x) 



dm 
dx 



x), 



and 



m 



X{x) dx 



(17.3) 



If A represents a function such as Ax'^ {0 < x < L), a bullwhip perhaps, then this is elementary, 
and TTitotai = I want to look at the reverse specification. Given an interval, I will specify the 

amount of mass in that interval and work backwards. The first example will be simple. The interval 
Xi < X < X2 iS denoted [a;i,X2]. The function m has this interval for its argument.* 



m{[xi,X2]) 



(17.4) 



(Xi < X2 < 0) 

(a;i < < a;2 < L) 
{xi<0<L< x-i) 
(0 < xi < X2 < L) 
(0 < xi < L < X2) 
[L <xi < X2) 

The density Ax"^ (0 < x < L) is of course a much easier way to describe the same distribution of mass. 





Ax'ys 

AL^/2, 
A{xl-x\)/2, 



A{L^ 




x?)/3 



Ax' 



on 



This distribution function, m([xi,X2]), comes from integrating the density function A(x 
the interval [xi,X2]. 

Another example is a variation on the same theme. It is slightly more involved, but still not too 

bad. 



m([xi,X2]) 



(Xi < X2 < 0) 

(xi < < X2 < L/2) 
(xi < < L/2 <X2<L) 
(xi < < L < X2) 
(0 < xi < X2 < L/2) 
(0 < xi < L/2 <X2<L) 
xf)/3 + mo (0 < xi < L/2 < L/2) 
(L/2 < xi < X2 < L) 
(L/2 < xi < L < X2) 
(L < Xi < X2) 





Axl/3 
Axl/?, + mo 
AL^/'i + mo 
A(x| -xf)/3 
A[x2 - x\)/3 + mo 
A(L3 
A{xl- 
A{L^ 




(17.5) 



xf)/3 
-x?)/3 



If you read through all these cases, you will see that the sole thing that I've added to the first example 
is a point mass mo at the point L/2. What density function A will produce this distribution? Answer: 
No function will do this. That's why the concept of a "generalized function" appeared. I could state 
this distribution function in words by saying 

"Take Eq. (17.4) and if [xi,X2] contains the point L/2 then add mo." 

That there's no density function A that will do this is inconvenient but not disastrous. When the very 
idea of a density was defined in Eq. (17.1), it started with the distribution function, the mass within 
the volume, and only arrived at the definition of a density by some manipulations. The density is a 
type of derivative and not all functions are differentiable. The function m([xi,X2]) or m(K) is more 
fundamental (if less convenient) than is the density function. 



* I'm abusing the notation here. In (17.2) m is a number. In (17.4) m is a function. You're used 
to this, and physicists do it all the time despite reproving glances from mathematicians. 
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17.2 Functionals 



/oo 
dxf{x)(j){x) 
-oo 



defines a scalar-valued function of a function variable. Given any (reasonable) function (f) as input, it 
returns a scalar. That is a functional. This one is a linear functional because it satisfies the equations 

F[a(P] = aF[(P] and + 02] = i^[0i] + i^[02] (17.6) 

This isn't a new idea, it's just a restatement of a familiar idea in another language. The mass density 
can define a useful functional {linear density for now). Given dm/dx = \{x) what is the total mass? 



/oo 
dx\{x)l = Mtotal 
-oo 



Where is the center of mass? 



1 



dx X{x)x = Xc 



^total J —oo 

Restated in the language of functionals, 

/■oo 

F[4>]= dxX{x)4){x) then M^otai = i^[l], ^cr. = ^ F[x] 

J-oo -'l^ total 

If, instead of mass density, you are describing the distribution of grades in a large class or the 
distribution of the speed of molecules in a gas, there are still other ways to use this sort of functional. 
If dN/dg = f{g) is the grade density in a class (number of students per grade interval), then with 

F[cP] = Jdgfig)cf>{g) 

A^students = F[l], mean grade = g = F[g], (17.7) 

-'''students 

variance = = F[{g - gf], skewness = Zs^K^ - 9?] 

-'''students ''students^' 

kurtosis excess = — — ^ 4^F[{g — g)'^] — 3 

-''students^ 

Unless you've studied some statistics, you will probably never have heard of skewness and kurtosis 
excess. They are ways to describe the shape of the density function, and for a Gaussian both these 
numbers are zero. If it's skewed to one side the skewness is non-zero. [Did I really say that?] The 
kurtosis excess compares the flatness of the density function to that of a Gaussian. 

The Maxwell-Boltzmann function describing the speeds of molecules in an ideal gas is at tem- 
perature T 

/MB(.) = (2^)'^'We— (178) 

dN/dv is the number of molecules per speed interval, but this function i^MB is normalized differently. 
It is instead {dN / dv) /Nx,ox.a\- It is the fraction of the molecules per speed interval. That saves carrying 
along a factor specifying the total number of molecules in the gas. I could have done the same thing in 
the preceding example of student grades, defining the functional Fi = .F/-/Vstudents- Then the equations 
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(17.7) would have a simpler appearance, such as g = Fi[g]. For the present case of molecular speeds, 
with F[(t)] = dv fMB{v)(t){v), 



m = 1, 



F[v] 



8kT 



nm 



F[mv^/2] = K.E. = ^kT 



;i7.9) 



Notice that the mean kinetic energy is not the kinetic energy that corresponds to the mean speed. 

Look back again to section 12.1 and you'll see not only a definition of "function" but a definition 
of "functional." It looks different from what I've been using here, but look again and you will see that 
when you view it in the proper light, that of chapter six, they are the same. Equations (12.3)-(12.5) 
involved vectors, but remember that when you look at them as elements of a vector space, functions 
are vectors too. 

Functional Derivatives 

In section 16.2, equations (16.6) through (16.10), you saw a development of the functional derivative. 
What does that do in this case? 



dx f{x)(f){x) 



so F[(f) + 6(f)] - F[(j)] = Jdx f{x)6(f){x) 



The functional derivative is the coefficient of 6(f) and dx, so it is 



(17.10) 



That means that the functional derivative of m(^[xi, X2\) in Eq. (17.4) is the linear mass density, 
X{x) = Ax"^, (0 < X < L). Is there such a thing as a functional integral? Yes, but not here, as it goes 
well beyond the scope of this chapter. Its development is central in quantum field theory. 

17.3 Generalization 

Given a function /, I can create a linear functional F using it as part of an integral. What sort of linear 
functional arises from /'? Integrate by parts to find out. Here I'm going to have to assume that / or 
(f) or both vanish at infinity, or the development here won't work. 



F[(/.]= / dxf{x)(^{x) 

J — oo 

dxf'{x)(l){x) = f{x)(j){x 



then 



dxf{x)(f)'ix) = -F[4>'] 



(17.11) 



In the same way, you can relate higher derivatives of / to the functional F. There's another restriction 
you need to make: For this functional — to make sense, the function (f) has to be differentiable. 
If you want higher derivatives of /, then needs to have still higher derivatives. 

If you know everything about F, what can you determine about /? 
If you assume that all the functions you're dealing with are smooth, having 
as many derivatives as you need, then the answer is simple: everything. If 
I have a rule by which to get a number F[(j)] for every (smooth) (j), then 
I can take a special case for and use it to find /. Use a that drops 
to zero very rapidly away from some given point; for example if n is large 
this function drops off rapidly away from the point Xq. 




(l>n{x) 
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I've arranged it so that the integral of (f)n over all x is one. If I want the value of / at Xq I can 
do an integral and take a limit. 

F[(j)n] = r dxf{x)<Pn{x)= r rfx/(x)./^e-"(^-^°)' 

J —oo J~oo V ^ 

As n increases to infinity, all that matters for / is its value at Xq. The integral becomes 

lim r rfx/(xo)A/^e-"(^-^o)' = f{xo) (17.12) 



This means that I can reconstruct the function / if I know everything about the functional F. To get 
the value of the derivative f'{xo) instead, simply use the function —0^ and take the same limit. This 
construction is another way to look at the functional derivative. The equations (17.10) and (17.12) 
produce the same answer. 

You say that this doesn't sound very practical? That it is an awfully difficult and roundabout 
way to do something? Not really. It goes back to the ideas of section 17.1. To define a density you 
have to know how much mass is contained in an arbitrarily specified volume. That's a functional. It 
didn't look like a functional there, but it is. You just have to rephrase it to see that it's the same thing. 

As before, do the case of linear mass density. A(x) is the density, and the functional F[(f)] = 
c^^; A(x)0(x). Then as in Eq. (17.4), m([xi,X2]) is that mass contained in the interval from Xi 
to X2 and you can in turn write it as a functional. 

Let y be the step function yix) = l \ /"^^i ~ ^. — /^^^ 
^ ^ ^ \0 (otherwise) 

/■X2 

then / dx X{x) = F[x\ = m{[xi,X2]) 

J X\ 

What happens if the function / is itself not differentiable? I can still define the original functional. 
Then I'll see what implications I can draw from the functional itself. The first and most important 
example of this is for / a step function. 

/(x) = ^(x) = I J F[<j)] = I^Jxe{x)<j){x) (17.13) 

9 has a step at the origin, so of course it's not differentiable there, but if it were possible in some way 
to define its derivative, then when you look at its related functional, it should give the answer —F[(f)'] 
as in Eq. (17.11), just as for any other function. What then is —F[(f)']7 

/oo /■oo 
dx9{x)(j)'{x) = - dx(j)'{x) = -^ix) =0(0) (17.14) 
-oo Jo 

This defines a perfectly respectable linear functional. Input the function and output the value of the 
function at zero. It easily satisfies Eq. (17.6), but there is no function / that when integrated against 
will yield this result. Still, it is so useful to be able to do these manipulations as if such a function 
exists that the notation of a "delta function" was invented. This is where the idea of a generalized 
function enters. 
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Green's functions 

In the discussion of the Green's function solution to a differential equation in section 4.6, I started with 
the differential equation 

mx + kx = F(t) 

and found a general solution in Eq. (4.34). This approach pictured the external force as a series of 
small impulses and added the results from each. 

'impulse 




x{t) = p dt'G{t-t')F{t') where ^(t) = | ^^^^o^ ^ -^^^ 

I wrote it a little differently there, but it's the same. Can you verify that this is a solution to the stated 
differential equation? Simply plug in, do a couple of derivatives and see what happens. One derivative 
at a time: 

^ = 1^ dt' Git - t')F{t') where G{t) = { '^'^^^ {J | (17.16) 

Now for the second derivative. Oops. I can't differentiate G. It has a step at t = 0. 

This looks like something you saw in a few paragraphs back, where there was a step in the 
function 6. I'm going to handle this difficulty now by something of a kludge. In the next section you'll 
see the notation that makes this manipulation easy and transparent. For now I will subtract and add 
the discontinuity from G by using the same step function 6. 

^[coscjot- 1 + 1] (t>0) 
(* < 0) 

Mcosu;ot-1] (J^O)Ul^(t) = Go(t) + -^(t) (17.17) 
(r < 0) J m m 



G{t) 



The (temporary) notation here is that Go is the part of G that doesn't have the discontinuity at t = 0. 
That part is differentiable. The expression for dx/dt now has two terms, one from the Gq and one 
from the 9. Do the first one: 

1^ dt' Goit - t')F{t') = dt' j^Goit - t')F{t') 

and ^GoW^{„^l-»^'-»'l (^12} 

The original differential equation involved mx + kx. The Go part of this is 

Hi' / ,ri[-^osina;o(t - t')] {t > t') \ 
\0 (t<f)J 

i-oo 10 {t<t')\ ^ > 



f oo 

m 

-oo 
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Use k = mujQ, and this is zero. 

Now go back to the extra term in 9. The kx terms doesn't have it, so all that's needed is 

mx + kx = / dt' —9{t - t')F{t') = ^ / dt' F{t') = F{t) 
dt m dt 

This verifies yet again that this Green's function solution works. 

17.4 Delta-function Notation 

Recognizing that it would be convenient to be able to differentiate non-differentiable functions, that it 
would make manipulations easier if you could talk about the density of a point mass (m/O =?), and 
that the impulsive force that appears in setting up Green's functions for the harmonic oscillator isn't a 
function, what do you do? If you're Dirac, you invent a notation that works. 

Two functionals are equal to each other if they give the same result for all arguments. Does that 
apply to the functions that go into defining them? If 

/oo 
dx f2ix)(f){x) 
-oo 

for all test functions (f) (smooth, infinitely differentiable, going to zero at infinity, whatever constraints 
we find expedient), does it mean that /i = /2? Well, no. Suppose that I change the value of /i at 
exactly one point, adding 1 there and calling the result f2- These functions aren't equal, but underneath 
an integral sign you can't tell the difference. In terms of the functionals they define, they are essentially 
equal: "equal in the sense of distributions." 

Extend this to the generalized functions 

The equation (17.14) leads to specifying the functional 

m = m (17.18) 

This delta-functional isn't a help in doing manipulations, so define the notation 

dx6{x)(f)ix) = 6[4)] = 0(0) (17.19) 



This notation isn't an integral in the sense of something like section 1.6, and 6{x) isn't a function, 
but the notation allows you effect manipulations just as if they were. Note; the symbol 5 here is not 
the same as the 5 in the functional derivative. We're just stuck with using the same symbol for two 
different things. Blame history and look at problem 17.10. 

You can treat the step function as differentiable, with 6' = 5, and this notation leads you 
smoothly to the right answer. 



Let 0^.„(x) = |j =^(^-^o), then j dxe^,{x)<P{x) = j 



dx (f){x) 



The derivative of this function is 



■^e^^^Xx) = s{x - Xo) 



You show this by 



oo 



dOxoix) 



dx — -J (f){x) = — dx6xo{x)(f)'{x) = — dx(j)'{x) = (f){xo) 

oo J —oo J Xo 
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The idea of a generalized function is that you can manipulate it as if it were an ordinary function 
provided that you put the end results of your manipulations under an integral. 

The manipulations for the harmonic oscillator in the previous section, translated to this language 
become 

sin Wot {t > 0) 
(t<0) 

Then the solution for a forcing function F{t) is 



mG + kG = 5{t) for G{t) 



muio 




xit) 



G{t - t')F{t') dt' 



because 



mx 



kx 



[m 



G + kG)F{t')dt' 



6{t - t')F{t') dt' = F{t) 



This is a lot simpler. Is it legal? Yes, though it took some serious mathematicians (Schwartz, Sobolev) 
some serious effort to develop the logical underpinnings for this subject. The result of their work is: It's 
o.k. 

17.5 Alternate Approach 

This delta-function method is so valuable that it's useful to examine it from more than one vantage. 
Here is a very different way to understand delta functions, one that avoids an explicit discussion of 
functionals. Picture a sequence of smooth functions that get narrower and taller as the parameter n 
gets bigger. Examples are 



IT ' 



n 



1 



1 sinnx 



n 



71 1 + n^X^ 71 X 



sech nx 



7T 



(17.20) 



Pick any one such sequence and call it 5n{x). (A "delta sequence") The factors in each case are 
arranged so that 

dx6n{x) = 1 



As n grows, each function closes in around x = and becomes very large there. Because these are 
perfectly smooth functions there's no question about integrating them. 



dx 6nix)4>{x) 



(17.21) 



makes sense as long as (j) doesn't cause trouble. You will typically have to assume that the (f) behave 
nicely at infinity, going to zero fast enough, and this is satisfied in the physics applications that we 
need. For large n any of these functions looks like a very narrow spike. If you multiply one of these 
6nS by a mass m, you have a linear mass density that is (for large n) concentrated near to a point: 
A(x) = m5n{x). Of course you can't take the limit as n — > oo because this doesn't have a limit. If 
you could, then that would be the density for a point mass: m5{x). 

What happens to (17.21) as n oo? For large n any of these delta- 
sequences approaches zero everywhere except at the origin. Near the origin 
(t){x) is very close to </)(0), and the function 5n is non-zero in only the tiny 
region around zero. If the function (p is simply continuous at the origin you 
have 




lim 

n— s>oo 



dx 6nix)(p{x) = 0(0) 



lim 

71— !>00 



dx 5nix) = (f)(0) (17.22) 
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At this point I can introduce a notation: 

dxS{x)(p{x)" MEANS 



lim 

71— >00 



dx 6nix)4>{x) 



(17.23) 



In this approach to distributions the collection of synnbols on the left has for its definition the collection 
of symbols on the right. Those in turn, such as J, have definitions that go back to the fundamentals of 
calculus. You cannot move the limit in this last integral under the integral sign. You can't interchange 
these limits because the limit of 5^ is not a function. 

In this development you say that the delta function is a notation, not for a function, but for a 
process (but then, so is the integral sign). That means that the underlying idea always goes back to a 
familiar, standard manipulation of ordinary calculus. If you have an equation involving such a function, 
say 



e'{x) = 6{x), 



then this means 



and that 



lim 



dx 9'^{x)(f){x) 



lim 



dx 6nix)(t>{x) = (f)(0) 



(17.24) 



-oo J —oo 

How can you tell if this equation is true? Remember now that in this interpretation these are sequences 
of ordinary, well-behaved functions, so you can do ordinary manipulations such as partial integration. 
The functions 6n are smooth and they rise from zero to one as x goes from — oo to +oo. As n becomes 
large, the interval over which this rise occurs will become narrower and narrower. In the end these On 
will approach the step function 9{x) of Eq. (17.13). 



dx 9n{x)(j){x) = 9n{x)(t>{x] 



dx 9n{x)(f)'{x) 



The functions go to zero at infinity — they're "test functions" — and that kills the boundary terms, 
leaving the last integral standing by itself. Take the limit as — )• oo on it. You can take the limit inside 
the integral now, because the limit of 9n is a perfectly good function, even if it is discontinuous. 



lim 

n 



dx9n{x)(j)' {x) 



dx9{x)(j)'{x) 



dx( 



(X 



X 



0(0) 



This is precisely what the second integral in Eq. (17.24) is. This is the proof that 9' = 6. Any proof of 
an equation involving such generalized functions requires you to integrate the equation against a test 
function and to determine if the resulting integral becomes an identity as n — > oo. This implies that 
it now makes sense to differentiate a discontinuous function — as long as you mean differentiation "in 
the sense of distributions." That's the jargon you encounter here. An equation such as 9' = 6 makes 
sense only when it is under an integral sign and is interpreted in the way that you just saw. 

In these manipulations, where did I use the particular form of the delta sequence? Never. A 
particular combination such as 



9n{x) 



1 r 2 

- 1 + - tan" 

2 L TT 



nx\ 



and 



5n{x) 



n 



1 



TT 1 + Tl^x"^ 



(17.25) 



never appeared. Any of the other sequences would have done just as well, and all that I needed was the 
properties of the sequence, not its particular representation. You can even use a delta sequence that 
doesn't look like any of the functions in Eq. (17.20). 



5n{x) 



2 /^gmx2 



(17.26) 
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This turns out to have all the properties that you need, though again, you don't have to invoke its 
explicit form. 

What is 5{ax)7 Integrate 5n{ax) with a test function. 

lim / dx5n{ax)(t){x) = Ymi I — 5n{y)(f){y/a) 

where y = ax. Actually, this isn't quite right. If a > it is fine, but if a is negative, then when 
X — )■ — oo you have y — )• +oo. You have to change the limits to put it in the standard form. You can 
carry out that case for yourself and verify that the expression covering both cases is 

f°° f°° dy 1 

lim/ dx Sniax)(j){x) = lim — 5n{y)(t){y / a) - 



\a\ 

Translate this into the language of delta functions and it is 

5{ax) = ^5(a;) (17.27) 
You can prove other relations in the same way. For example 

5(x2-a2) = ^[5(x-a) + 5(x + a)] or 5{f{x))=Y,^jr^\^i^-^k) (17.28) 

In the latter equation, X}^ is a root of /, and you sum over all roots. Notice that it doesn't make any 
sense if you have a double root. Just try to see what 5{x'^) would mean. The last of these identities 
contains the others as special cases. Eq. (17.27) implies that 5 is even. 

17.6 DifTerential Equations 

Where do you use these delta functions? Start with differential equations. I'll pick one that has the 
smallest number of technical questions associated with it. I want to solve the equation 

^-k'f = F{x) (17.29) 

subject to conditions that f{x) should approach zero for large magnitude x. I'll assume that the given 
function F has this property too. 

But first: Let k and y be constants, then solve 



1-2 

dx"^ 



k^g = 5{x-y) (17.30) 



I want a solution that is well-behaved at infinity, staying finite there. This equality is "in the sense of 
distributions" recall, and I'll derive it a couple of different ways, here and in the next section. 

First way: Treat the 6 as simply a spike at x = y. Everywhere else on the x-axis it is zero. Solve 
the equation for two cases then, x < y and x > y. In both cases the form is the same. 

g" - k^g = 0, so g{x) = Ae^"" + Be-^"" 

For X <y, \ want g{x) to go to zero far away, so that requires the coefficient of e~^^ to be zero. For 
X > 0, the reverse is true and only the e^^^ can be present. 

( \ — \ "^^^^ < y) 

~ \ Be~k^ {x > y) 
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Now I have to make g satisfy Eq. (17.30) at x = y. 

Compute dg/dx. But wait. This is impossible unless the function is at least continuous. If it 
isn't then I'd be differentiating a step function and I don't want to do that (at least not yet). That is 

g{y-) = Ae^y = g{y+) = Be-^v (17.31) 
This is one equation in the two unknowns A and B. Now differentiate. 



(17.32) 



cf^ ^ f Ake^^ (x < y) 
dx \ —Bke~^^ {x > y) 

This is in turn differentiable everywhere except at x = y. There it has a step 

discontinuity in g' = g'{y+) - g'{y-) = -Bke'^^ - Ake^v 

This means that (17.32) is the sum of two things, one differentiable, and the other a step, a multiple 
of 0. 

^ = differentiable stuff + ( - Bke'^y - Ake^y)e{x - y) 
The differentiable stuff satisfies the differential equation g" — k'^g = 0. For the rest. Compute d'^g/dx^ 

( - Bke-'^y - Ake^y) ^e{x -y) = {- Bke-^y - Ake^y)S{x - y) 

Put this together and remember the equation you're solving, Eq. (17.30). 

g" - k'^g = (from the differentiable stuff) 

+ {-Bke-^y -Ake^y)6{x-y) =5{x~y) 

Now there are two equations for A and B, this one and Eq. (17.31). 

Ae'^y = Be'^y A = -e'^y /2k 

solve these to get 

-Bke-^y - Ake^y = 1 B = -e^y/2k 

Finally, back to g. 



When you get a fairly simple form of solution such as this, you have to see if you could have saved 
some work, perhaps replacing brute labor with insight? Of course. The original differential equation 
(17.30) is symmetric around the point y. It's plausible to look for a solution that behaves the same 
way, using {x — y) as the variable instead of x. Either that or you could do the special case y = and 
then change variables at the end to move the delta function over to y. See problem 17.11. 
There is standard notation for this function; it is a Green's function. 




17 — Densities and Distributions 



12 



Now to solve the original differential equation, Eq. (17.29). Substitute into this equation the 
function 



dyG{x,y)F{y) 



dyG{x,y)F{y) 

d^Gix,y) 



dy 



dx"^ 



PG{x,y) 



F{y) 



/oo 
dy5[x-y)F{y)=F{x) 
-oo 



This is the whole point of delta functions. They make this sort of manipulation as easy as dealing with 
an ordinary function. Easier, once you're used to them. 

For example, if = Fq between —xq and +Xq and zero elsewhere, then the solution for / is 



/ 



dyG{x,y)F{y) 



dyFQe-^\''-y\/2k 



Fo 
2k 

Fo 
fc2 



' J^l^dye-^^'^-y^ {x>xo) 

I-^jco e"^^^"^^ + dy e-'^^y-'^^ {-xq < x < xq) 
J%dye-''^y-'') [x<-xq) 



g-fca; gjj^]^ ^2;q (x > Xq) 

[1 — e"^^" coshfca;] (— xq < a; < xq) 
e^^ sinh kxn 



(17.35) 



'^0 



{X < —Xq) 



You can see that this resulting solution is an even function of x, necessarily so because the original 
differential equation is even in x, the function F{x) is even in x, and the boundary conditions are even 
in X. 

Other Difrerential Equations 

Can you apply this method to other equations? Yes, many. Try the simplest first order equation: 



dG 

dx 
d[_ 
dx 



5{x) 
9{x) 



G{x) = e{x) 



dx'G{x — x')g{x') 



dx' g{x') 



which clearly satisfies df /dx = g. 

If you try d'^G/dx'^ = S{x) you explain the origin of problem 1.48. 

Take the same equation d'^G/dx'^ = 6{x — x') but in the domain < x < L and with the 
boundary conditions G'(O) = = G{L). The result is 



G{x) 



x{x'~L)/L (0 < X < x') 
x'{x-L)/L {x' < X < L) 




(17.36) 



17.7 Using Fourier Transforms 

Solve Eq. (17.30) another way. Fourier transform everything in sight. 

/oo 
dx Six - y)e-'i'' (17.37) 
-oo 



d^ 

dx"^ 



k'^g = 5{x-y) 



dx 



d^ 
dx"^ 



k\ 
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The right side is designed to be easy. For the left side, do some partial integrations. I'm looking for a 
solution that goes to zero at infinity, so the boundary terms will vanish. See Eq. (15.12). 



/oo 
dx[- q"^ - k^]g{x)e~^'^^ = e" 
-oo 

The left side involves only the Fourier transform of g. Call it g. 

[-q^-k']giq) = e-'iy, so ~g{q) -- 
Now invert the transform. 



iqy 



(17.38) 



-iqy 



q^ + k"^ 



g{x) 



dq 
2tx 



iqx 



dq e^i^^-y) 
27C q^ + k'^ 



Do this by contour integration, where the integrand has singularities at g = ±.ik. 



Ci 
— yj- 



1 r Qiq{x-y) 



The poles are at ±.ik, and the exponential dominates the behavior of the integrand at large \q\, so 
there are two cases: x > y and x < y. Pick the first of these, then the integrand vanishes rapidly as 

q — > +ioo. 



'Co 




(5^ 



C4 



Res 



2n q=ik A;2 + g2 



Compute the residue. 



Jq{x-y) 



Jq{x-y) 



Jq{x-y) 



k'^ + q^ {q — ik){q + ik) {q — ik){2ik) 
The coefficient of l/(g — ik) is the residue, so 

^-k{x-y) 



2k 



[x > y) 



(17.39) 



in agreement with Eq. (17.33). The x < y case is yours. 



17.8 More Dimensions 

How do you handle problems in three dimensions? 5{f) = 5{x)5{y)5{z). For example I can describe 
the charge density of a point charge, dq/dV as q5{r — Tq). The integral of this is 



q6{f— fo) d^r = q 
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as long as the position fo is inside the volume of integration. For an example that uses this, look at 
the potential of a point charge, satisfying Poisson's equation. 

VV = -p/eo 

What is the potential for a specified charge density? Start by finding the Green's function. 

y'^G = 5{r), or instead do: V'^G - k'^G = 5{r) (17.40) 

The reason for starting with the latter equation is that you run into some problems with the first form. 
It's too singular. I'll solve the second one and then take the limit as A; — > 0. The Fourier transform 
method is simpler, so use that. 

j (Pr [V'G - k^G] ■ ^ = 1 

When you integrate by parts (twice) along each of the three integration directions dx, dy, and dz, you 
pull down a factor of —q^ = — g| — Qy — qI just as in one dimension. 

j dh[- q^G - k^G] e'"^- ^=1 or G{q) = -^^^ 



where G is, as before, the Fourier transform of G. 

Now invert the transform. Each dimension picks up a factor of l/27r. 



d^q e 



iq • r 



|3 g2 _|_ /j2 



This is a three dimensional integral, and the coordinate system to choose is spherical, doesn't 
depend on the direction of q, and the single place that this direction matters is in the dot product in 
the exponent. The coordinates are q, 6, (j) for the vector q, and since I can pick my coordinate system 
any way I want, I will pick the coordinate axis along the direction of the vector r. Remember that in 
this integral f is just some constant. 

"1^ f- giqr cos 9 

G = - / q'^dq sin 9 d9d(f) 



The integral d(f) becomes a factor 27r. Let u = cos6, and that integral becomes easy. All that is left is 
the dq integral. 



^_ r tdq 1 r iqr 



More contour integrals: There are poles in the g-plane at g = ±ik. The q^ factors are even. 
The q in the denominator would make the integrand odd except that the combination of exponentials 
in brackets are also odd. The integrand as a whole is even. I will then extend the limits to ±oo and 
divide by 2. 



Q = ^ [°° ^ fg*?^ _ 
SnHr y_oo q'^ + k^^ 
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There are two terms; start with the first one. r > 0, so e*'^''" as g — +ioo. The contour is along 
the real axis, so push it toward too and pick up the residue at +ik. 

' ^ - ' :2vr^ Res , ,,. e--^- = -J^Jmi^e-'^- (17.41) 



Svr^ir J Srr'^ir q=ik {q — ik){q + ik) 871'^ir 2ik 

For the second term, see problem 17.15. Combine it with the preceding result to get 

G = ^e"^^ (17.42) 

Anr 

This is the solution to Eq. (17.40), and if I now let k go to zero, it is G = — l/47rr. 

Just as in one dimension, once you have the Green's function, you can write down the solution 
to the original equation. 

VV = -p/eo =^ V = ~- [ d^r'G{f,f')p{f') = [ d^r' (17.43) 

Co J 47reo J \r -r'\ 

State this equation in English, and it says that the potential of a single point charge is q/Aneor 
and that the total potential is the sum over the contributions of all the charges. Of course this 
development also provides the Green's function for the more complicated equation (17.42). 

Applications to Potentials 

Let a charge density be qS{f). This satisfies / d'^r p = q. The potential it generates is that of a point 
charge at the origin. This just reproduces the Green's function. 

, d\' q6(f')—^ = (17 AA) 

47reo J \r -r'\ Aneor 

What if p(r ) = —pd6{f)/dz7 [Why —p7 Patience.] At least you can say that the dimensions 
of the constant p are charge times length. Now use the Green's function to compute the potential it 
generates. 



4,71 eo J dz' \r — f'\ Aireo J dz'\f—r'\ 

p d 1 



r '=0 



47160 dz' \ f—r' 

This is awkward, so I'll use a little trick. Make a change of variables in (17.45) u = r — f' , then 

d d p d I 



(17.45) 



p d I 
f =0 4^^o duz u 



dz' duz ' 47reo dz' \r — r 

(See also problem 17.19.) Cutting through the notation, this last expression is just 

. ~p d 1 —p d 1 p z 



(17.46) 



Aireo dz r 47reo dz -\-y'^ + z'^ 47reo (3,2 + _^ ^2)3/2 

p z p cos 6 



(17.47) 
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The expression —diXjr^jdz is such a simple way to compute this, the potential of an electric 
dipole, that it is worth trying to understand why it works. And in the process, why is this an electric 
dipole? The charge density, the source of the potential is a derivative, and a derivative is (the limit of) 
a difference quotient. This density is just that of two charges, and they produce potentials just as in 
Eq. (17.44). There are two point charges, with delta-function densities, one at the origin the other at 
-lA^. 



[5(r + 5 Az) — (5(r )] gives potential 



-V 



47reoAz 



1 



1 



Jr + f Az| 



(17.48) 





The picture of the potential that arises from this pair of charges is (a). A negative charge (— g = 
— p/ A^) at —z^z and a corresponding positive charge +g at the origin. This picture explains why this 
charge density, p(f ) = —pd6{f)/dz, represent an electric dipole. Specifically it represents a dipole 
whose vector representation points toward +2, from the negative charge to the positive one. It even 
explains why the result (17.47) has the sign that it does: The potential (f) is positive along the positive 
2;-axis and negative below. That's because a point on the positive ^-axis is closer to the positive charge 
than it is to the negative charge. 

Now the problem is to understand why the potential of this dipole ends up with a result as simple 
as the derivative {—p/4:neQ)d{l/r) /dz: The potential at the point P in the figure (a) comes from the 
two charges at the two distances ri and r2. Construct figure (b) by moving the line r2 upward so that 
it starts at the origin. Imagine a new charge configuration consisting of only one charge, q = p/ Az at 
the origin. Now evaluate the potentials that this single charge produces at two different points Pi and 
P2 that are a distance Az apart. Then subtract them. 



</'(P2)-0(Pl 



In the notation of the preceding equations this is 



0(P2)-0(Pl) 



47reo 



1 

r2 



1 

ri 



47reo [|^+ zAz\ 



P 



AneoAz 



\r + zAz\ r 



Except for a factor of (—1), this is Eq. (17.48), and it explains why the potential caused by an ideal 
electric dipole at the origin can be found by taking the potential of a point charge and differentiating 
it. 

q , d 



-'point charge 



dz 



-'point charge 



47reor 

Here, the electric dipole strength is qa. 

Can you repeat this process? Yes. Instead of two opposite point charges near each other, you 
can place two opposite point dipoles near each other. 



^linear quadrupole 



d_ 
dz 



(17.49) 
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This is the potential from the charge density p = +Q d'^6{r)/dz^ , where Q = qa?. [What about 
d"^ /dxdyl] 

Exercises 

1 What is the analog of Eq. (17.4) for the linear mass density \{x) = C (a constant) for < x < L 
and zero otherwise? 

2 Take the preceding mass density and add a point mass mo at a; = L/2. What is the distribution 
m{[xi,X2]) now? 

3 Use the A from the first exercise and define the functional F[<j)\ = J^^dx X{x)(j){x). What is the 
total mass, F[l] = Ml What is the mean position of the mass, F[x\/Ml 

4 As in the preceding exercise, what are the variance, the skewness, and the kurtosis excess? 

5 What is dx 5{x — Xq)1 

6 Pick any two of Eq. (17.20) and show that they are valid delta sequences. 

7 What is f^^dt6{t)7 

8 In Eq. (17.12) the function (j)n{x) appeared. Sketch a graph of it and of — 0'^(x), which is needed 
in Eq. (17.11). 
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Problems 

17.1 Calculate the mean, the variance, the skewness, and the kurtosis excess for a Gaussian: f{g) = 
j\^Q-B{g-go) (_oo < g < oo). Assume that this function is normalized the same way that Eq. (17.8) 
is, so that its integral is one. 

17.2 Calculate the mean, variance, skewness, and the kurtosis excess for a flat distribution, f{g) = con- 
stant, (0 < 5f < S-max)- Ans: Var = gl,/l2 kurt. exc. = -6/5 

17.3 Derive the results stated in Eq. (17.9). Compare mv'^ /2 to K.E. Compare this to the results of 
problem 2.48. 

17.4 Show that you can rewrite Eq. (17.16) as an integral J^^dt' cos Uo{t — t')F(t') and differen- 
tiate this directly, showing yet again that (17.15) satisfies the differential equation. 

17.5 What are the units of a delta function? 

17.6 Show that 

5{f{x))=5{x-xo)/\f'ixo)\ 

where Xq is the root of /. Assume just one root for now, and the extension to many roots will turn 
this into a sum as in Eq. (17.28). 

17.7 Show that 

(a) x6'{x) = -6{x) (b) x5(x) = 

(c) d'{-x) = -6'{x) (d) f{x)6{x -a) = f{a)6{x - a) 

17.8 Verify that the functions in Eq. (17.20) satisfy the requirements for a delta sequence. Are they 
normalized to have an integral of one? Sketch each. Sketch Eq. (17.26). It is complex, so sketch both 
parts. How can a delta sequence be complex? Verify that the imaginary part of this function doesn't 
contribute. 

17.9 What is the analog of Eq. (17.25) if 6n is a sequence of Gaussians: ^Jn / ir e^'^^^ 1 
Ans: 9nix) = 5 [l + erf [xy/n]] 

17.10 Interpret the functional derivative of the functional in Eq. (17.18): S 5[(j)]/S(j). Despite appear- 
ances, this actually makes sense. Ans: 5{x) 

17.11 Repeat the derivation of Eq. (17.33) but with less labor, selecting the form of the function g to 
simplify the work. In the discussion following this equation, reread the comments on this subject. 

17.12 Verify the derivation of Eq. (17.35). Also examine this solution for the cases that Xq is very 
large and that it is very small. 

17.13 Fill in the steps in section 17.7 leading to the Green's function for g" — k'^g = 5. 

17.14 Derive the analog of Eq. (17.39) for the case x < y. 

17.15 Calculate the contribution of the second exponential factor leading to Eq. (17.41). 
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17.16 Starting with the formulation in Eq. (17.23), what is the result of 5' and of 5" on a test function? 
Draw sketches of a typical Sn, 5'^, and 5'^. 

17.17 If p(r) = qa? d'^5{r)/dz^, compute the potential and sketch the charge density. You should 
express your answer in spherical coordinates as well as rectangular, perhaps commenting on the nature 
of the results and relating it to functions you have encountered before. You can do this calculation in 
either rectangular or spherical coordinates. 

Ans: {2qa?/A'KtQ)P2{cose)/r^ 

17.18 What is a picture of the charge density p{f) = qa?' d'^5{r) / dxdyl 
(Planar quadrupole) What is the potential for this case? 

17.19 In Eq. (17.46) I was not at all explicit about which variables are kept constant in each partial 
derivative. Sort this out for both d/dz' and for d/duz- 

17.20 Use the results of the problem 17.16, showing graphs of 5n and its derivatives. Look again at 
the statements leading up to Eq. (17.31), that g is continuous, and ask what would happen if it is not. 
Think of the right hand side of Eq. (17.30) as a 5n too in this case, and draw a graph of the left side 
of the same equation if Qn is assumed to change very fast, approaching a discontinuous function as 
n — > oo. Demonstrate by looking at the graphs of the left and right side of the equation that this can't 
be a solution and so that g must be continuous as claimed. 

17.21 Calculate the mean, variance, skewness, and the kurtosis excess for the density f{g) = A[5{g) + 
5{g — go) + 6{g — xgo)]. See how these results vary with the parameter x. 

Ans: skewness = 2^^/2(1 + x)(x - 2)(2x - l)/(l - x + x^) 
kurt. excess = -3 + f (l + x*^ + (1 - x)^)/(l - x + x'^)^ 

VJ .22 Calculate the potential of a linear quadrupole as in Eq. (17.49). Also, what is the potential of 
the planar array mentioned there? You should be able to express the first of these in terms of familiar 
objects. 

17.23 (If this seems out of place, it's used in the next problems.) The unit square, < x < 1 and 
< ?/ < 1, has area f dxdy = 1 over the limits of x and y. Now change the variables to 

u=^{x + y) and v = x — y 

and evaluate the integral, [ dudv over the square, showing that you get the same answer. You have 
only to work out all the limits. Draw a picture. This is a special example of how to change multiple 
variables of integration. The single variable integral generalizes 



/f dx f [ u\X v) 

f{x)dx = J f{x)-^du to j f{x,y)dxdy = J /(x, y) ^^^' ^^ dudv 



where 



(9(x, y) _ / dx/du dx/dv 
d{u, v) ~ ^ \dy/du dy/dv 



For the given change from x, y to u,v show that this Jacobian determinant is one. A discussion of the 
Jacobian appears in many advanced calculus texts. 

17.24 Problem 17.1 asked for the mean and variance for a Gaussian, f{g) = Ae"^^^"^"-*^ . Interpreting 
this as a distribution of grades in a class, what is the resulting distribution of the average of any two 
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students? That is, given this function for all students, what is the resulting distribution of {gi +g2)/27 
What is the mean of this and what is the root-mean-square deviation from the mean? How do these 
compare to the original distribution? To do this, note that f{g)dg is the fraction of students in the 
interval g to g + dg, so f{gi)f{g2)dgi dg2 is the fraction for both. Now make the change of variables 

X =^(91+92) and y = g^-g2 
then the fraction of these coordinates between x and x + dx and y and y + dy \s 



f{9i)f{92)dxdy = f{x + y/2)f{x - y/2)dxdy 

Note where the result of the preceding problem is used here. For fixed x, integrate over all y in order 
to give you the fraction between x and x + dx. That is the distribution function for {gi + 5'2)/2. 
[Complete the square.] Ans: Another Gaussian with the same mean and with rms deviation from the 
mean decreased by a factor \/2. 



17.25 Same problem as the preceding one, but the initial function is 

= i-oo<g< 00) 

In this case however, you don't have to evaluate the mean and the rms deviation. Show why not. 
Ans: The result reproduces the original f{g) exactly, with no change in the spread. These two problems 
illustrate examples of "stable distributions," for which the distribution of the average of two variables 
has the same form as the original distribution, changing at most the widths. There are an infinite 
number of other stable distributions, but there are precisely three that have simple and explicit forms. 
These examples show two of them. The Residue Theorem helps here. 

17.26 Same problem as the preceding two, but the initial function is 

(a) /(^) = for 0<^<^max (b) /(^) = J - ^1 ) + - ^2)] 



17.27 In the same way as defined in Eq. (17.10), what is the functional derivative of Eq. (17.5)? 

17.28 Rederive Eq. (17.27) by choosing an explicit delta sequence, 5n{x). 

17.29 Verify the result in Eq. (17.36). 



